Four-coloring model and frustrated superfluidity in the diamond lattice 
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We propose a novel four-coloring model which describes "frustrated superfluidity" of p-band 
bosons in the diamond optical lattice. The superfluid phases of the condensate wavefunction on 
the diamond-lattice bonds are mapped to four distinct colors with the constraint that four differ- 
ently colored bonds meet at the same site. We show that the combined effects of strong intra-site 
repulsions, inter-site phase coherence, and anisotropy give rise to a macroscopically degenerate 
ground state at the classical level. By further mapping the problem to an antiferromagnetic Potts 
model on the pyrochlore lattice, we demonstrate that both color and orbital angular momentum 
correlations exhibit power-law decay in the degenerate manifold that is described by an emergent 
magnetostatic theory with three independent flux fields. 

PACS numbers: 03.75.Ss, 05.50.+q, 71.10.Fd, 73.43.Nq 



Strongly frustrated systems are hosts to various com- 
plex orders, unusual phases and elementary excitations. 
A well-studied example is the emergence of a critical 
Coulomb phase in systems ranging from water ice [l[, 
pyrochlore magnets (spin ice) [2|, [3| to p-band fermions 
in the optical diamond lattice (orbital ice) [4]. A ubiq- 
uitous feature here is the appearance of an extensively 
degenerate manifold due to a large number of uncon- 
strained degrees of freedom. Although the macroscopic 
degeneracy is usually of geometrical origin, especially for 
magnetic systems [5|, the nontrivial spatial dependence 
of anisotropic orbital interactions in optical lattices pro- 
vides a new playground for studying interesting phenom- 
ena related to a highly degenerate manifold [4j, l6|4S| . In 
particular, for p-orbital bosonic condensates in optical 
lattices, the frustrated couplings between the U(l) phase 
degrees of freedom can lead to intriguing quantum many- 
body states P, [To|| . 

For conventional superfluid states of bosons, their 
ground state wavefunctions are positive-definite as stated 
by the "no-node" theorem which is valid under very 
general conditions It means that the superfluid 

phases are uniform over the entire condensate, in other 
words, non-frustrated. Exotic states of bosons beyond 
the "no-node" theorem have been theoretically proposed 
and experimentally realized in ultra-cold atom optical 
lattices [12] . For example, when bosons are pumped 



to high orbital bands, their condensate wavefunctions 
form non-trivial representations of the lattice symmetry 
group, dubbed "unconventional" Bose-Einstein conden- 
sations, say, the complex- valued p + ip-type condensates 
with spontaneously broken time-reversal symmetry (see 
Ref. [9] for a brief review). These condensates are meta- 
stable excited states, and hence are not constrained by 
the "no-node" theorem. 

In this article we study the Bose-Einstein condensates 
formed by bosons pumped into the p-orbital bands in 
a diamond lattice. We show that the problem of inter- 



site phase coherence in this lattice represents a highly 
frustrated system and can be mapped to a four-coloring 
model on the diamond lattice. By combining analytical 
arguments and numerical simulations, we demonstrate 
the existence of an algebraic Coulomb phase in the highly 
degenerate ground state of the diamond-lattice 4-coloring 
model. Our work thus provides a three-dimensional (3D) 
generalization of the celebrated three-coloring model on 
the planar honeycomb lattice [T^j, and a similar 4- 



coloring problem on the square lattice [14| . In both cases 
the ground state is extensively degenerate and exhibits 
power-law color correlations. 

Model Hamiltonian. We begin with a discussion of the 
Hamiltonian for p-orbital bosons in the diamond optical 
lattice. This bipartite lattice can be generated by the in- 
terference of four laser beams with suitably arranged light 
polarizations as discussed in Ref. [15]. The unit vectors 
from each site in sublattice A to its four neighboring B- 
sites are denoted as no = [111], ni = [111], 112 = [111], 
and 113 = [111]. Around the center of each lattice site, 
the point group symmetry is of which the degenerate 
Px,y,z~OTbitals form a triplet irreducible representation. 
The band energy is thus represented as 



v=0 <y>||n„ 



(1) 



where p v is the projection of the p-orbitals along the 
direction of h u : p v = p x h% + p y h v v + p z h z v . We only 
keep the longitudinal hopping between p-orbitals along 
the bond direction, i.e., the a-bonding, and neglect the 
small transverse hopping terms. Here tu is positive be- 
cause of the odd parity of p-orbitals. 

The general on-site interaction has the form: 

#int = T,i V ab;cd pt,a P{b Pi,cPi,d> where V a b;cd = 

g J d 3 r ?/j a (r)^(r)^ c (r)'0d(r), g is the contact interac- 
tion parameter, and ip a (a = x, y, z) are Wannier func- 
tions of p-orbitals. The Td group contains the symmetry 
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FIG. 1: (Color online) Mapping between the orbital angular 
momentum and color configuration on the diamond lattice. 
The red, green, blue, and yellow bonds have phases = 0°, 
90°, 180°, and 270°, respectively, when projected to the base 
plane perpendicular the direction of angular momentum L. 
The mapping is not one-to-one, as each L corresponds to 4 
different cyclic permutations of the colors which preserve the 
same chirality. 

of rotation around the z-axis followed by the inversion, 
which transforms ip z — >• —i/j z ,iJj x — » —ip y , and ip y — >• ip x . 
Since V a b ;c d should be invariant under this transforma- 
tion, symmetry consideration indicates that there are 
only two independent parameters V xx;xx and V xy;xy . The 
interaction terms can thus be reorganized as 

H int = Uj2{nl-l^}+AUY: £ n? )Q , (2) 

i i a=x, y ,z 

where is the particle number in p a -orbital (a = 
x,y,z) at site r^, rti = rii jX + n^y + n^ z , and = 
— \tyiv\V%V\ is the orbital angular momemtum opera- 
tor. The first term in Eq. (|2]) represents the Hubbard 
interaction with the spherical symmetry around the site 
center and U = 3V xy , xy as given in Ref. @, E3 ■ We 
assume repulsive interactions U > 0. The ferro-orbital 
interaction means that bosons prefer to maximize the 
on-site orbital angular momentum such that their wave- 
functions are most extended spatially to reduce the onsite 
repulsion [9j, analogous to the second Hund's rule of elec- 
trons filling in degenerate atomic orbit als. The second 
term with AU = \(V XX - XX — SV xy ^ xy ) comes from the 
symmetry crystal field. Compared with the spherically 
symmetric case (AU = 0), the angular distributions of 
Wannier orbitals ip a in the diamond lattice expand from 
±x,±y?iz toward the the bond directions h u . This en- 
hances the inter-orbital repulsion but weakens the intra- 
orbital one, such that AU < 0. 

Mapping to the 4~ coloring model We consider the 
strong coupling limit with a weak crystal field, i.e., 
U ^> ty and U ^> \AU\. Assuming that the average filling 
is large, we approximate each site as a small condensate 



and treat it classically. In order to minimize the U -term, 
bosons condense into the orbital of [p^(ei) -\- i ($2)] / V% 
with ei _L e2 at each site, where p^(e^) is the projection of 
the p-orbitals on the direction of e$. The corresponding 
angular momentum is L || ei x &2. 

The crystal field from a AU < term aligns L to di- 
rections along ±x, ±y, =tz. Assuming the average filling 
per site is n and neglecting the quantum depleting of the 
condensation fraction, the kinetic energy of each nearest- 
neighbor bond becomes 

2 

Eij = -^ nt \\ cos (</v ~ <A?>) » ( 3 ) 

where (ij) \\ h v , <p iiV is the superfluid phase along bond 
(ij) at site r^, and the factor | comes from the angles 
between the bond direction h u and the base planes of 
ei 5 2 for Li at site i and Lj at site j, respectively. The 
U(l) phase in each site winds 2tt around the nodal line 
which is parallel to L^. 

The phase 0^ along the bond direction h u is the same 
as that along the projection of h u in the ei^ plane; see 
Fig. [TJ For AU < 0, the ei^ plane is perpendicular to 
one of the cubic axes. Since the projections of the four 
h u in the xy, yz and zx planes are evenly distributed, 
(j)i yU take values of up to an overall U(l) phase, 

where — 0, 1, 2, 3 defines the color indices for the bond 
(ij) || h u . They correspond to colors R, G, B, and Y, 
respectively, in Fig. [TJ We can thus express <pi^ as 

7T 

<t>i i v=0i-\-Gij—, (4) 

where 6i denotes the U(l) phase fluctuation. The con- 
straint of phase coherence is that each bond can only be 
assigned with one color and four bonds connecting each 
site take different colors (Fig. [T]). 

To minimize the kinetic energy, we freeze 0i = for 
all sites, such that ^ v = (f)j yU for each bond. The classic 
ground states of the p-band bosons thus map into a highly 
frustrated 4-coloring model. An example of such config- 
urations in the diamond lattice is shown in Fig. [2J There 
are 4 ! = 24 coloring configurations on each site z, and 6 
different orientations (±x, d=y, d=z) for L^, which are de- 
fined as "chirality" of the condensate. For each chirality, 
there are 24/6 = 4 coloring configurations corresponding 
to a C4 rotation around the direction of L^; see Fig. [TJ 
Consequently, the classical ground states correspond to 
those of the 4-coloring model modular 4. The chirality 
can be explicitly computed as || (n* — no) x (n^ — n^), 
where denotes the bond with phase tp at site r^. 

The phase correlation of the p-band condensate has 
contributions from both the discrete color variables err- 
and the U(l) variables 0i. At the classic level, it is dom- 
inated by the color correlations at low temperatures in 
which the gapped excitations of defects can be neglected. 
For all the defect-free configurations, the partition func- 
tions for thermal phase fluctuations are identical, and 
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FIG. 2: (Color online) Coloring of the dimaond lattice that 
satisfy the constraints that four bonds attached to the same 
vertex have different colors. By assigning a 4- state Potts vari- 
able to label the coloring of the bond, the 4-coloring model on 
the diamond lattice can be mapped to an antiferromagnetic 
Potts model on its medial pyrochlore lattice. 
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FIG. 3: (Color online) (a) The angular averaged color correla- 
tion function C(r) on a system with linear size L — 10 in the 
ground state of the pyrochlore 4-state Potts model. The inset 
shows the four unit vectors r v (y — ~ 3) corresponding to 
the four different colors, (b) Entropy S and specific heat C as 
a function of temperature T obtained from the Monte Carlo 
simulations. The energy is measured in units of J 7 . The en- 
tropy S(T) is obtained by numerically integrating the specific 
heat curve; the integration constant So = 0.2112. 



thus the color correlations of and phase fluctuations of 
0i are decoupled. In 3D below a critical temperature T c , 
there exists a long-range-order for the U(l) phase 0$, thus 
up to a temperature-dependent factor, the correlations of 
(f)i yU are just described by the color correlations. However, 
quantum phase fluctuations may select an optimal color- 
ing pattern through the "order by disorder" mechanism, 
which should occur at a very low temperature scale T' c . 
Thus the 4-coloring model applies at « T « T c , and 
we will see that algebraic phase correlation develops in 
this temperature regime. 

Antiferromagnetic Potts model on pyrochlore lattice. 
Since the color variables are defined on bonds of the dia- 
mond lattice, it is more convenient to consider an equiv- 
alent Potts model in the pyrochlore lattice whose sites 
correspond to the bond midpoints in the diamond lattice 
(Fig. [2}. A Potts variable <j m = aij = 0, 1, 2, 3 is assigned 
to each pyrochlore site such that the corresponding phase 
is given by Eq. (j4]). The constraint that bonds attached 
to the same vertex have different colors translates to an 
antiferromagnetic interaction between nearest-neighbor 
Potts variables on the pyrochlore: 



H 



j ^2 ^ ar > 

(ml) 



(5) 



Here the coupling constant J ~ nt\\. To avoid con- 
fusion, we use m, I to label the pyrochlore sites. The 
hard constraint of the 4-coloring model is recovered in 
the J — >• oo limit. Consider first a single tetrahedron, 
there are 4 ! = 24 ground states in which the four sites 
are in different Potts states. These correspond to the 
24 coloring schemes on a diamond site. When general- 
izing to the pyrochlore lattice, model (j5j) represents an 
under-constrained system with an extensively degenerate 
ground state, similar to their two-dimensional counter- 
parts [II EI]. 



In order to determine the structure of the degenerate 
manifold, we employ Monte Carlo simulations with non- 
local loop updates on a finite lattice. In each loop update, 
two sites are chosen randomly in a given ground state. 
These two sites will necessarily be of different color, say 
R and B. With periodic boundary conditions, a RB- 
colored loop containing the two chosen sites is uniquely 
determined. By exchanging colors R and B over the 
length of the loop, the non-local update results in a new 
ground state as all the color constraints remain satis- 
fied. There are a total of L) = 6 types of loops with 
different colors. Since all the microstates have equal sta- 
tistical weight, loop updates of all lengths and colors are 
accepted in order to satisfy detailed balance. 

We first study the color correlations in the ground 
state. To this end, we introduce four unit vectors t v 
pointing toward different corners of a regular tetrahedron 
as shown in the inset of Fig. [3](a). The correlation func- 



tion between two sites is defined as C(r m i) = (r 



which equals 1 for a m = a\ and —1/3 for a m ^ g\. 
Fig. [3]( a) shows the color correlation averaged over differ- 
ent angular directions obtained from Monte Carlo simula- 
tions. The correlations C(r) fall off rapidly beyond a few 
nearest-neighbor distances, indicating a color-disordered 
ground state. 

The degeneracy of the ground state can be computed 
by a method similar to the Pauling estimation for the 
degeneracy of water ice [l|. Consider a single tetrahe- 
dron, 4 ! out of 4 4 Potts configurations satisfy the color 
constraint. Treating the constraints imposed by differ- 
ent tetrahedra as independent, the number of degenerate 
ground states is W ~ 4 Ar (4!/4 4 ) Ar / 2 , where N is the num- 
ber of pyrochlore sites and N/2 is the number of tetrahe- 
dra (or diamond sites). This estimation gives an entropy: 
S /Nk B = (l/2)ln(3/2) « 0.2027. The residual entropy 
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can also be estimated with the aid of Monte Carlo sim- 
ulations. The system entropy as a function of the tem- 
perature can be obtained by numerically integrating the 
specific- heat curve. The integration constant So at zero 
temperature can be fixed by the fact that the high-T en- 
tropy density of the Potts model should approach In 4. 
As shown in Fig. [3(b), the numerically determined resid- 
ual entropy So/Nks ~ 0.2112 agrees very well with the 
Pauling estimate. 

Emergent Coulomb phase. Although the Potts ground 
state is disordered, the local hard constraints give rise to 
nontrivial long-range correlations, similar to the case in 
dimer models |17| and spin ice 18, ljj. To make manifest 
the analogy with frustrated spin systems, we consider a 
Heisenberg antiferromagnet in the pyrochlore lattice: 



(6) 



Here J' — 3 J/4 and S m denotes a classical 0(3) vector of 
unit length. The second term represents a special 'tetra- 
hedral' anisotropy: the function /(S) has four degenerate 
minimum at the tetrahedral vectors t v . In the D — >> oo 
limit, Hamiltonian © reduces to the Potts model as the 
spins are aligned to the tetrahedral vectors according to 
S m = r am . Moreover, since the exchange term in (j6]) can 
be rewritten as (J'/2)^2^ |S^| 2 , where denotes the 
total spin of a tetrahedron, a Potts configuration that 
satisfies the 4-color constraint thus maps to a ground 
state of the spin model as = X^=o Tv ~ 0- 

The mapping to Heisenberg spins also allows us to re- 
cast the color constraint into a conservation law of effec- 
tive magnetic flux [18-20]. We first define three 'mag- 
netic' fields, each corresponds to a component of the 
Heisenberg spin, at the diamond sites: 



u=0 



v=0 



(7) 



Here a = x,y,z and the pyrochlore site m corresponds 
to the bond (ij) || h u in the diamond lattice, i.e. 
&ij = S m . In the coarse-grained approximation, the color 
constraints = translate to a divergence constraint 
V • B a (r) = for the magnetic fields. 

The effective free energy of the ground-state mani- 
fold arises entirely from entropy and has the form of the 
magnetostatic theory with three independent flux fields 
F[B a (r)] oc E a /d 3 r|B a (r)| 2 . Essentially, it states 
that the partition function is dominated by microstates 
characterized by B a w 0. This is due to a large num- 
ber of flippable loops in such states. Although superfi- 
cially F describes Gaussian fluctuations of the magnetic 
fields, the divergence- free constraint in momentum space 
k • B a (k) = indicates that only transverse fluctuations 
are allowed. Consequently, the correct correlators are 
obtained by projecting out the longitudinal fluctuations. 




FIG. 4: (Color online) (a) Correlation function vs linear sys- 
tem size L from Monte Carlo simulations. Here C(L) de- 
notes the color correlation between two points separated by 
(L/2, L/2, 0) in a finite lattice containing L 3 cubic unit cells, 
(b) The length distribution P(s) of flippable loops from Monte 
Carlo simulations on a L = 12 system. 



The asymptotic field correlators in real-space has the fa- 
mous dipolar form 



<B?(r)Bf(0)> ex 5 aP {5 ab - 3f f b )/r 3 



(8) 



characteristic of a Coulomb phase [21 1. 

We also performed Monte Carlo simulations to inves- 
tigate the color correlators at large distances. Fig. [3(a) 
shows the correlation function C(L) between two sites 
separated by r = |f (1,1,0) as a function of the linear 
system size L. A power-law dependence C(L) ~ L~ 3 at 
large L confirms the emergence of the Coulomb phase 
in the Potts model. Given its critical nature, flippable 
loops (with alternating colors) of all sizes are expected 
to exist in the Coulomb phase. In numerical simulations, 
however, the largest loop is bounded by the system size. 
Fig. [3(b) shows the length distribution P(s) of loops from 
Monte Carlo simulations. Similar to those found in the 



Coulomb phase of spin ice 22], the distribution exhibits 
two regimes with distinct behaviors. For short loops, the 
probability function scales as P(s) ~ s -3 / 2 , while an s- 
independent regime is observed for larger s. 

In terms of the flux fields, the orbital angular momen- 
tum can be expressed as L = B a x B^, where the two 
components a ^ /3 depend on the particular mapping 
between phases and colors. The correlation function of 
angular momentum thus exhibits a power-law 1/r 6 decay 
in the degenerate manifold. 

Mapping for AU > 0. At last, if we take AU as a 
free variable and consider the situation of AU > 0, the 
classic ground states are also described by the 4-coloring 
model but with a different mapping. For each site i, 
the particle density in the three p-orbitals are equal in 
order to minimize the term of AC/, thus its angular mo- 
mentum hi is parallel or anti-parallel to one of the four 
bond directions. Such a bond is the nodal line of the 
onsite condensate, thus it is broken as marked by Y. 
The kinetic energy of any of the other three bonds (ij) 
becomes = — -nty cos^j, — (j)j yU ) where (ij) \\ h u . 
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Again <\>^ v can be decomposed as 0^ = 6^ + |cr^-7r, where 
dij = 0,1,2 are the color indices (marked as R, G, B, 
respectively). The even and odd permutations of the 3 
colors determines whether || h u or — h u . The 24 col- 
oring patterns for the four bonds connecting to site i is 
classified as 8 different chiralities, and for each chirality 
there are 3 configurations corresponding to a Cs rotation 
around L^. 

Conclusion. We have studied a diamond-lattice four- 
coloring model which describes the frustrated couplings 
between the superfluid phase degrees of freedom for p- 
band Bose-Einstein condensates in the diamond lattice. 
We have also shown that the ground states of the 4- 
coloring model is macroscopically degenerate and is de- 
scribed by an effective magnetostatics theory with three 
independent flux fields. Both color and orbital angular 
momentum correlations decay algebraically in this emer- 
gent Coulomb phase. Interestingly, point defects violat- 
ing the color constraints carry "magnetic" charges corre- 
sponding to two of the three flux fields. A future direc- 
tion of study is to explore the kinematics and dynamics 
of these novel quasiparticles. Finally, another challeng- 
ing task is to examine how the huge degeneracy is lifted 
through quantum order-by-disorder mechanism. 
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Note added. Upon completion of this work, we became 
aware of similar results on the pyrochlore 4-state Potts 
model in Ref. 0. 
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